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Low-energy electron impact excitations of N2 molecules are studied using the fixed-bond R-matrix 
method based on state-averaged complete active space SCF orbitals. Thirteen target electronic states of 
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N2 are included in the model within a valence configuration interaction representations of the target states. 
Integrated as well as differential cross sections of the A 3 X+, B 3 U g , W 3 A U , B a' a}H g , w 1 A u and 
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C 3 n„ states are calculated and compared with the previous experimental measurements. These excitations, 
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especially of the higher four states, have not been studied enough theoretically in the previous literature. 
In general, good agreements are observed both in the integrated and differential cross sections. However, 
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I. INTRODUCTION 



Electron impact excitation of nitrogen molecules plays an important role in atmospheric emis- 
sion of planets and satellites such as the Earth, Titan and Triton. For example, excitation of the 
a l Tlg state and subsequent transitions to the ground X 1 !^ state are responsible for the far ultravio- 
let emissions of the Lyman-Birge-Hopfield system which is prominent in the airglow of the Earth's 
atmosphere! 1]. Recently, Khakoo et al.J^I measured differential cross sections (DCSs) of electron 
impact excitation of N2 molecule from the ground X 1 !** state to the 8 lowest excited electronic 
states of A 3 Z+, B 3 Il g , W 3 A U , 5 ,3 £~, a' 1 !,', a l Ii g , w 1 A u and C 3 II l( states. Based on their differen- 
tial cross section data, Johnson et al.[3] derived integral cross sections (ICSs) for these electron 
impact excitations. In general, their ICSs are smaller than the other experimental cross sections at 
low impact energies below 30 eV. These deviations may have some significance on study of atmo- 
spheric emissions, because a mean kinetic energy of electron at high altitudes is about 10 eV|4|]. 
To shed light on this situation from a theoretical point of view, we perform the ab initio R-matrix 
calculations of electron impact excitations of N 2 molecule in this work. 

Many previous experimental measurements have been focused on excitation to a specific elec- 
tronic state. For example, Ajello and Shemansky[5] and Mason and Newell|^] measured ICSs for 



electron impact excitation to the a l U g state, whereas Poparic et 



al.fl, 



Zubek[8] and Zubek and 



King||9|] measured cross sections for the C 3 n„ state. In addition to these works, Zetner and Trajmar 
Jioll reported excitation cross sections to the A 3 E^, B 3 U g , W 3 A U and a l Ii g states. So far, compre- 
hensive measurements of the excitation to the 8 lowest electronic states are limited to three groups 

n n 

of Cartwright et al. [11], Brunger and Teubner [12] and Khakoo et al.[|2|]. The measurements of 
Brunger and TeubnerJl21 include excitation DCSs for the E 3 !^ and a" l ^Z states in addition to 
the 8 lowest excited states. The DCSs of Brunger and TeubnejJlJ] and Khakoo et al.]2] were later 
converted to ICSs by Campbell et al.[13] and Johnson et al.yfl, respectively. Detailed reviews on 
electron N 2 collisions can be found in Itikawa and Brunger and Buckman lll5n . 

Several groups have performed theoretical calculation of low energy electron collisions with 
N 2 molecule. For example, Chung and Lin[16] employed the Born approximation to calculate 
excitation cross sections for the 11 target states including the A 3 ?,*, B 3 U g , W 3 A U , a l Ii g , w 1 A u and 



C 3 n„ states. Later, the same group of Holley 



et al.Jl7n 



calculated excitation ICSs for the a l Tl<, state 



using a two-state-close-coupling method. Fliflet et al. lll8l1 and Mu-Tao and McKoyl I9I1 reported 
distorted-wave cross sections for excitation of the A 3 E+, B 3 U g , W 3 A U , w l A u , C 3 n„, E 3 2, g , b' T£ 
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and c' states. In general, these approximate methods are expected to be accurate at high impact 
energies above 30 eV. However, more elaborate method is required for precise comparison with 
experiment at low energies. Gillan et al.[20] calculated excitation ICSs for the A 3 Z^, B 3 Tl g and 
W 3 A U states using the fixed nuclei R-matrix method. They included the 4 lowest target states in 
their R-matrix model, with target CI wave functions containing 2-13 CSFs. Their cross sections 
for the A 3 Z+ and W 3 A U states agree well with the experimental results of Cartwright et al. Ill ill . 
However, ICSs for the B 3 U g state deviate considerably from the experimental cross sections. Sub- 



sequently, they extended their R-matrix model to include the 8 lowest valence states H2 111 . Their 
target CI wave functions were much improved from their previous work by employing valence 
active space description, resulting in 68-120 CSFs per target state. In their paper, the ICSs were 
shown for the A 3 2+, B 3 Tl g , W 3 A U and B ,3 "L~ states, while the DCSs were presented for only the 
A 3 £+ state. Agreement with the ICSs of Cartwright et al. ill ill is good for these 4 excited states. 
However, agreement is marginal at DCS level. 

In this work, we study electron impact excitation of N 2 molecule by the fixed nuclei R-matrix 
method as in our previous work on electron O2 scatterings 11221. l23ll . Although theoretical treat- 
ment is similar to the previous work of Gillan et al. Clll . more target states and partial waves of a 
scattering electron are included in the present work. Main purpose of this work is comparison of 
ICSs as well as DCSs for the 8 lowest excited states with the experimental results of Cartwright 
et al. Jill . Brunger and Teubner[12], Campbell et al.[13], Khakoo et al. [2] and Johnson et al.isl. 
This is because previous theoretical works have covered only a part of these 8 excitations. 

In this paper, details of the calculation are presented in section 2, and we discuss the results in 
section 3 comparing our ICSs and DCSs with the previous theoretical and available experimental 
data. Then summary is given in section 4. 
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II. THEORETICAL METHODS 



3,H 



as well 



The R-matrix method itself has been described extensively in the literature H24I . 
as in our previous paper !122ll . Thus we do not repeat general explanation of the method here. 
We used a modified version of the polyatomic programs in the UK molecular R-matrix codes [24]. 
These programs utilize the gaussian type orbitals (GTO) to represent target electronic states as well 
as a scattering electron. Although most of the previous R-matrix works in electron N 2 collisions 
have employed Slater type orbitals (STO), we select GTO mainly because of simplicity of the input 
and availability of basis functions. In the R-matrix calculations, we have included 13 target states; 
X 1 ^, A 3 I H + , B 3 U g , W 3 A ln B'%, a'%, a l Il g , wM H , C 3 U U , £ 3 Z;, a"^ g , c l U u and c' 1 ^. The 
potential energy curves of these target electronic states are shown in figure Q] for reference. These 
target states were represented by valence configuration interaction wave functions constructed by 
state averaged complete active space SCF (SA-CASSCF) orbitals. Note that some target states, 
E 3 ?,^, a" and c' 1 ?.*, are Rydberg states and cannot be described adequately in the present 
valence active space. Inclusion of these states are intended to improve quality of the R-matrix 



calculations by adding more targ et states in the model, as in our previous works 11221 . 



as other R-matrix works [27, 



230 as well 



2811 . Test calculation was performed with an extra Aa g orbital in the 



target orbital set. However, the target excitation energies as well as the excitation cross sections 
did not change much compared to the results with valence orbital set described above. Also, 
removal of 3b\ u orbital from target active space did not affect the result much in our calculation. 
In this study, the SA-CASSCF orbitals were obtained by calculations with MOLPRO suites of 
programs Ii29ll . The target orbitals were constructed from the [5s3pld] level of basis set taken 
from Sarpal et al. I130ll . In our fixed-bond R-matrix calculations, the target states were evaluated 
at the equilibrium bond length R = 2.068 a of the N 2 X 1 !,^ ground electronic state. Although we 
also performed calculations with R = 2.100 ao as in the previous R-matrix calculation of Gillan et 
al. [|2in . the cross sections with R = 2.068 a and R = 2.100 a are almost the same. Thus, we will 
only show the results with the equilibrium bond length of N 2 in the next section. The radius of the 
R-matrix sphere was chosen to be 10 a in our calculations. In order to represent the scattering 
electron, we included diffuse gaussian functions up to / = 5, with 9 functions for I = 0,1 functions 
for /= 1 - 3 and 6 functions for 1 = 4 and 5. Exponents of these diffuse gaussians were fitted using 
the GTOBAS program 113 ill in the UK R-matrix codes. In addition to these continuum orbitals, we 
included 8 extra virtual orbitals, one for each symmetry. 
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We constructed the 15-electron configurations from the orbitals listed in table HI The CI target 
wave functions are composed of the valence orbitals in table U with the la g and Ibu, orbitals kept 
doubly occupied. We have included 3 types of configurations in the calculation. The first type of 
configurations has the form, 

\al\b\ u {2a g 3a g lb 2 Ab, u 2b lu 3b u Ah g \b 2g ) w (1) 

here we assume that the total symmetry of this 15 electrons system is 2 A g . The first 4 electrons 
are always kept in the la g and \b\ u orbitals, then the next 10 electrons are distributed over the 
valence orbitals with restriction of target state symmetry, l A g symmetry of the N 2 ground state 
in this case. The last electron, the scattering electron, occupies one of the diffuse orbitals, of a g 
symmetry in this example. To complete the wave function with the total symmetry 2 A g , we also 
have to include configurations with the other target states combined with diffuse orbitals having 
appropriate symmetry in the same way as in the example. The second type of configurations has 
the form, 

l4l^J2apa g l& 2H l^ H 2& lH 3^ lH l^l^} 10 ( 1 A g ){4a,} 1 ( 2 A g ), (2) 

where the scattering electron occupies a bound 4a g extra virtual orbital, instead of the diffuse 
continuum orbitals in the expression (1). As in table U we included one extra virtual orbital for 
each symmetry. The third type of configurations has the form, 

la z g lbi u {2a g 3a g lb 2u lb 3u 2b lu 3b lu lb 3g lb 2g } n ( 2 A g ) . (3) 

In this case, the last 1 1 electrons including the scattering electron are distributed over the valence 
orbitals with the restriction of 2 A g symmetry. Note that the third type of configurations are crucial 
in description of N 2 resonance states, which often have dominant contributions to the excitation 
cross sections. In this way, the number of configurations generated for a specific total symmetry 
is typically about 60000, though the final dimension of the inner region Hamiltonian is reduced to 



be about 600 by using CI target contraction and prototype CI expansion method [32]. 

The R-matrix calculations were performed for all 8 irreducible representations of the D 2 /, sym- 
metry, A g , B 2u , B 3u , Bi g , Biu, B 3g , B 2g and A u , in doublet spin multiplicity of the electron plus target 
system. DCSs were evaluated in the same way as in our previous paper l!23n . 
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HI. RESULTS AND DISCUSSION 



A. Excitation energies 



Figure Q] shows the potential energy curves of all N 2 target states included in the present R- 
matrix model. These curves were obtained by the same SA-CASSCF method employed in our 
R-matrix calculation. Table |n] compares the excitation energies of the N 2 target states from the 
present calculation with the previous R-matrix results of Gillan et al.l 2lj . multi-reference coupled 
cluster results of Ben-Shlomo and Kaldor fl3 3M as well as experimental values. Since these ener- 
gies are evaluated at different inter-nuclear distance, 2.068 a in our case, 2.100 a in Gillan et 
al. ll2lll and 2.074 a in Ben-Shlomo and Kaldor 113 311 . precise comparison is not so meaningful. 
However, deviations of excitation energies from the experimental values are less than 0.8 eV in 
our calculation, which is good considering the level of calculation. In terms of excitation energies, 
our calculation and the previous R-matrix calculation of Gillan et al.[21] have similar quality. 

In addition to this good agreement of target energies with experimental results, energies 
are also well described in our SA-CASSCF calculation. In our calculation, X 2 ?.* and A 2 II l( 
states are located at 15.63 and 17.21 eV above N 2 X 1 !,^ state, respectively. Compared to the 
experimental values of 15.61 and 17.08 eV, our SA-CASSCF calculation gives good results. Note 
that the energy ordering of NJ X 2 ?.* and A 2 U U states are not well described in the Hartree Fock 
level calculation, see Ermler and McLean [34] for example. 



B. Integral cross sections 



Figure [2] shows integral cross sections for electron impact excitation from the N 2 X 1 ^ state 
to the A 3 E+, B 3 U g , W 3 A U and B' 3 1,~ states. In this figure, present results are compared with the 
previous R-matrix calculations of Gillan et al. 112 ill , recent calculations of da Costa and Lima [35], 



experimental results of Cartwright et al. (|ll|] . Campbell et alJLj] and recent measurements of 



Johnson et al. Dll. Renormalized values of Cartwright et al. ill ill are used as recommended by 
Trajmar et al. 113611 . Figure [3] compares the present excitation cross sections of the a' x YT u , a l Tl g , 



w l A u and C 3 II i( states with the previous experimental results of Cartwright et al. ill ill . Campbell 



et al.[13] and Jo 



inson et al.[3]. For the a'FL state cross sections, the recent calculations of da 



Costa and Lima [35], other experimental values of Ajello and Shemansky [5], Zetner and Trajmar 
Jlpj and Mason and Newell [6] are included. For the C 3 n„ state cross sections, the experimental 



results of Zubek [8|], Zubek and King [Q] and Poparic et al. [7] are included. 

Our excitation cross sections for the A 3 £+ state have a resonance feature at approximately 
12 eV as in the previous R-matrix results of Gillan et al. 112 ill . The 2 n„ resonance state is 
responsible for this peak structure. The main configuration of this resonance state is ln 3 u ln 2 . 
Other than the 2 H U symmetry partial cross sections, The 2 Ii g symmetry contributes to the ICSs 
as a smooth background component (not shown in the figure). Compared to the previous R- 
matrix cross sections, the peak at 12 eV is more pronounced in our case. Our results are slightly 
larger than theirs at 12-17.5 eV. Compared to the recent experimental results of Johnson et aLyD, 
our cross sections are about 50% larger at 12.5-20 eV, though 50% smaller at 10 eV. Also our 



calculation overestimates the results of Campbell et al. lll3h . however, the results of Cartwright et 
al. lllin agree well with our results except at 12.5 eV. The position of the resonance peak depends 
rather strongly on the inter-nuclear distance of N 2 molecule, which is 12.2 eV for 2.068 Oq and 
1 1.75 eV for 2.100 a in our calculations. Thus, inclusion of vibrational motion may be necessary 
to resolve this discrepancy of the resonance peak. 

Our excitation cross sections for the B 3 U g state have a small bump at 12.8 eV, which is not 
evident in the previous R-matrix cross sections. The origin of this bump is the l 2 A g state, with 
main configuration of 3cr g \n 2 . Other than this bump, the ICSs are mostly composed of the 2 Ii g 
symmetry contribution and have a shape similar to the previous R-matrix results. The magnitude of 
our ICSs is about 50% larger than the previous results of Gillan et al.[21]. Recently, da Costa and 
Lima [1350 calculated ICSs for the B 3 Xl g state using the Schwinger multichannel method with the 
minimal orbital basis for the single configuration interactions (MOB-SCI) approach. There cross 
sections are much larger than our results above 12 eV. Also, there is a prominent peak around 10 
eV in their ICSs, which does not exist in the R-matrix calculations. Compared to the experimental 
ICSs, our results agree well with the cross sections of Cartwright et al. [11], especially above 15 
eV. However, the results of Campbell et al. ll 1311 are much larger than ours. Recent measurements 



of Johnson et al.yO agree better with the previous R-matrix calculation of Gillan et al.pll]. 

For the excitation cross sections for the W 3 A U state, our results have a shape and magnitude 
similar to the previous R-matrix results. Most of our ICSs are composed of the 2 U g symmetry 
partial cross sections. Agreement with the experimental cross sections of Johnson et al. fl3|] is good 
in this case. The cross sections of Campbell et al. II 1311 agree well with our results at 15 and 17.5 
eV, but their value is about half as much as our result at 20 eV. The results of Cartwright et al. Jlli 
are about two times larger than our cross sections. 
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Our excitation cross sections for the B' 3 T.~ state are about half of the previous R-matrix cross 
sections of Gillan et al.[21]. Apart from this difference in magnitude, the shape of the cross 
sections is similar. Dominant component in these ICSs is the 2 U g symmetry partial cross sections, 
although the 2 U U symmetry also has certain contribution around 18-20 eV. Among 3 different 



experimental measurements, our results agree well with the results of Johnson et al.[3]. The 
experimental cross sections of other two groups are much larger than our results at 15 and 17.5 eV, 
and have a different energy dependence compared to the present calculation. 

The situation of the excitation cross sections for the a' 1 E~ state is similar to the case of the 
5' 3 2~ state. The 2 H g and 2 JT„ symmetry partial cross sections contribute almost equally to the 
ICSs. Our cross sections roughly agree with the results of Johnson et al.yQ, while the cross 
sections of Cartwright et al. lilD and Campbell et al. [|13|l at 15 eV are much larger than our result. 
The results of Cartwright et al. [11] and Campbell et al.||13fl decrease as impact energy increases 
from 15 to 20 eV, however, our cross sections increase mildly in this energy region. 

In case of excitation to the a 1 Tl g state, several other experimental results are available in addition 



3fl, and Johnson et al.Ol . The cross 



to the measurements of Cartwright et al. [11], Campbell et al.[ 
section profiles of Johnson et al.yfl, Ajello and Shemansky [|5|], Cartwright et al. Hi] an d Mason 
and Newell [6] are similar to our ICSs. However, the magnitude of our cross sections is lower 
than the experimental values in most case except the cross sections of Johnson et al.yO. At 15, 
17.5 and 20 eV, agreement of our results with the cross sections of Johnson et al.[3] is very good, 
although our cross section at 12.5 eV is twice as large as their value. Note that there is no dominant 
symmetry contribution to the calculated ICSs. All partial cross sections contribute rather equally 
to the ICSs. Recent ICSs of da Costa and Lima [35] by the Schwinger multichannel method are 
also shown in the panel (b) of the figure [3] Their result has a sharp peak at 12 eV as in their 
calculation for the B 3 U g state excitation. This difference between our and their results may come 
from different number of target states considered in the scattering calculation. Only the X l l. g , 
a 1 IIg and B 3 U g states were included in the calculations of da Costa and Lima. The other part of 
the cross section profile is similar to the shape of our cross sections, although the magnitude of 
their cross sections are about twice as large as our results at 15-20 eV. 

Our excitation cross section for the w l A u state gradually increases as a function of energy 
from the threshold to the broad peak around 17.5 eV, then decreases toward 20 eV. In this case, 
agreement with the results of Johnson et al. Q3|] is not so good compared to the excitations of the 
a x Ti g and a' l L~ states. Our cross sections are about 50% larger than their values at 17.5 and 20 eV. 
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At 15 eV, our results agree well with the cross section of Johnson et al-Jsj], however, they are about 
50% lower than the results of Cartwright et al. [11] and Campbell et al. II 1311 . In the calculated 
ICSs, the 2 n„ symmetry partial cross section is a major component, with a minor contribution 
from the 2 Ii g symmetry. 

The calculated excitation cross sections for the C 3 n„ state has a peak similar to the experimental 
results of Zubek [8] and Poparic et al.Q7[|. Although the shape of the cross sections is similar, 
position of the cross section peak is different from experimental results. In our case, it is located at 
about 17 eV, whereas corresponding peaks are located at 14 eV in the experimental cross sections. 
The height of the peak in our ICSs is lower than the experimental values of Zubek [8] and Poparic 
et al.[7]. It is unclear whether there is a cross section peak in the experimental cross sections of 
Cartwright et al. Ill ill . Campbell et al. 113(1 and Johnson et al.[3]. At least, it appears that they do 
not have a peak around 17 eV. The origin of this discrepancy in the cross section peak is uncertain, 
but may be related to the employment of the fixed-nuclei approximation or insufficiency of higher 
excited target states in the R-matrix model. The calculated ICSs are composed of the 2 E+ and 
2 £~ symmetry partial cross sections near the peak structure at 17 eV. The contribution of the 2 £+ 
symmetry is about 50% larger than the 2 Z~ component. Other than these two symmetries, the 2 n g 
symmetry partial cross section contributes to the ICSs as a smooth background component. 



C. Differential cross sections 



Figure H] shows calculated DCSs for excitation of the A 3 2+ state with the experimental results 
of Khakoo et al. Brunger and Teubner 1 121. Cartwright et al. [ 11], Zetner and Trajmar Jioll . 



LeClair and Trajmar [37] and the previous R-matrix DCSs of Gillan et al. 112 ill . Our DCSs at 12.5, 
15 and 17.5 eV have similar shape in common. They are enhanced in backward direction and 
have a small dimple at 120 degrees with a bump at 75 degrees. At 17.5 eV, our cross sections 
are located between the experimental values of Khakoo et al.[2] and Cartwright et al. |ll|] . The 



profile of the experimental DCS s are reproduced well in our calculation. At 1 5 eV, our results agree 
better with the results of Khakoo et al. [2] compared to the other experiments. In the DCSs of 



the previous R-matrix calculation of Gillan et al.[21], a bump is located at 40 degrees and a small 
dimple is located at 100 degrees, which agree better with the experimental results of Brunger and 
Teubner[12]. In our calculation, these dimple and bump are shifted toward backward direction by 
20 degrees, and agreement with the results of Brunger and Teubner [12] is not so good. At 12.5 
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eV, our calculation overestimates the experimental results by a factor of two. As seen in panel (a) 
of Figf2l this discrepancy is related to the existence of a resonance peak around 12.5 eV. 

Figure [5] compares calculated excitation DCSs for the B 3 U g state with the experimental and 
recent theoretical results. Our DCSs at 12.5, 15 and 17.5 eV have backward-enhanced feature 
with a broad peak at 130 degrees. At 15 and 17.5 eV, our DCSs agree well with the results of 
Khakoo et al. yO at forward direction below 80 degrees. However, their DCSs are smaller than 
ours by a factor of two at 80-130 degrees. Agreement with the results of Cartwright et al. Ill IN at 
15 eV is good at 20- 130 degrees, although their DCSs are twice as large as our DCSs at 17.5 eV for 
low scattering angles. Because of a resonance-like feature at 12.5 eV as seen in panel (b) of Fig[2l 
our results are larger than the experimental results at 12.5 eV. Recent Schwinger multi-channel 
results of da Costa and Lima [35] are much larger than our DCSs at 12.5 and 15 eV. The deviation 
is especially large at 12.5 eV, which is possibly related to the difference in the excitation energies 
of the target state. 

Figure [6] shows the excitation DCSs for the W 3 A U state with the experimental cross sections. 
At 15 and 17.5 eV, our cross section gradually increases as a function of scattering angle, without 
noticeable bump or dip. At 12.5 eV, the shape of DCSs is nearly symmetric around 90 degrees. 
Agreement with the experimental DCSs of Khakoo et al. [2] is good, although their results at 
15 and 17.5 eV have more complex structure such as a small peak at 80 degrees. Our DCSs are 
generally smaller than the other experimental results of Brunger and Teubner lil2ll . Cartwright et 
al. Ill], Zetner and Trajmar | lp| . 

Excitation cross sections for the B' 3 "L~ state are shown in figure [TJ Calculated DCSs decrease 
to be zero toward and 180 degrees, because of a selection rule associated with transition 



L38, 



3911 . Our DCSs have a broad single peak near 90 degrees at 12.5 and 15 eV, whereas there are 



two broad peaks at 17.5 eV. The position of the right peak at 17.5 eV coincides with that of the 
experimental DCSs of Khakoo et al. |2J] and Cartwright et al. [11], although the peak of Cartwright 
et al. Ill ill is much higher than ours. Our results agree well with the DCSs of Khakoo et al. Q2Q at 
15 and 17.5 eV. However, their cross sections at 15 eV have a small dip at 100 degrees and a small 
bump 60 degrees, which do not exist in our results. At 12.5 eV, our cross sections are slightly 
larger than the results of Khakoo et al.Q2[|. On the whole, agreement with the other experimental 
results of Brunger and Teubner Ill2n and Cartwright et al. Ill ill is not good. 

Figure [8] shows the excitation DCSs for the a' x YT u state. Because of selection rule, DCSs 

at and 180 degrees become zero as in the case of the B' 3 1,~ state DCSs. Calculated DCSs have a 
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broad single peak near 60 degrees at 12.5 and 15 eV. At 17.5 eV, there are two broad peaks at 50 
and 120 degrees. Although there is slight overestimation of DCSs near 50-60 degrees, our DCSs 
agree marginally with the results of Khakoo et al.|20- Agreement with the other experimental 
results is not good except low scattering angles at 17.5 eV. 

Figure |9] compares our excitation DCSs for the a l U g state with the experimental cross sec- 
tions. Because of large variation of the DCSs, the cross sections are shown in logarithmic scale. 
Calculated DCSs are strongly forward-enhanced, which is consistent with all experimental results 
shown in the figure. Our DCSs at 12.5 eV have a small dip around 100 degrees, which moves 
forward to 85 degrees at 15 eV and 75 degrees at 17.5 eV. This behavior roughly agrees with the 
results of Cartwright et al ii ill and Khakoo et al.yt]. At 15 eV, our DCSs agree better with the 
results of Khakoo et al.[Q] than the other experimental DCSs. At 17.5 eV, the results of Cartwright 
et al.[ll] are closer to our DCSs at scattering angles above 40 degrees. Below 40 degrees, our 
calculation significantly underestimates the experimental DCSs. Our results at 12.5 eV are located 
between the DCSs of Cartwright et al. jll|l and Khakoo et al.Q2Q, however the shape of the DCSs 
is similar to their results. The shapes of DCSs calculated by da Costa and Lima 13511 are similar to 
our results. However, their cross sections are larger than our results at low- scattering angles below 
80 degrees, where their results agree better with the experimental DCSs of Brunger and Teubner 
lll2ll and Zetner and Trajmar [10]. 

Figure [10] shows calculated excitation DCSs for the w l A u state with the experimental cross 
sections. Our DCSs are enhanced in forward direction as in the case of the a l U g state. However, 
magnitude of the enhancement is much smaller than that of the a 1 ^ state. Agreement with the 
DCSs of Cartwright et al. \lu\ is good at 17.5 eV except low scattering angles below 20 degrees. 
At 12.5 and 15 eV, their results are much larger than our DCSs. At 15 eV, our DCSs agree 
marginally with the results of Khakoo et al. [2], although details of the DCS profile are different. 
Their results are smaller than ours at 17.5 and 12.5 eV. Discrepancy is especially large for forward 
scattering at 12.5 eV. 

Figure \TT\ shows excitation DCSs for the C 3 n„ state with the experimental cross sections of 
Khakoo et al.[2], Brunger and Teubner [|12|l . Zubek and King [9] and Cartwright et al. Cal- 
culated DCS profiles are almost fiat at 12.5 and 15 eV, whereas they are enhanced in backward 
direction at 17.5 eV. Below 90 degrees, slope of the calculated DCSs at 17.5 eV is similar to the 
results of Khakoo et al.]2], Zubek and King [9] and Cartwright et al. Jlli . though our results are 
about 50% larger than their DCSs. In general, our results do not agree well with the experimental 
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DCSs. Although the ICS of Khakoo et al. [2] at 15 eV agrees well our result as shown in panel (d) 
of Figf3l angular dependence of the cross sections appears to be different. 



D. Discussion 



The excitation ICSs of the B 3 U g state, shown in panel (b) of figure |2l have a small bump around 



13 eV. However, there is no such structure in the previous R-matrix ICSs of Gillan et al. 112 ill . The 
origin of this bump in our calculation is the l 2 A g state, with main configuration of 3<j g \n 2 . 
The existence of the N~ l 2 A g state can also be verified by usual CASSCF calculation on N~ with 
valence active space ignoring continuum orbitals. In molpro calculation, the energy of the l 2 A g 
state is 15.7 eV. Since diffuse continuum orbitals are added in the R-matrix calculation, the energy 
of the state is stabilized to be 12.8 eV in the present scattering calculation. In the same way, the 
2 H U (Itt^Itt 2 ,) resonance peak in the A 3 2+ excitation ICSs can be verified by usual CASSCF 
calculation. In molpro calculation, it is located at 14.7 eV, whereas the position of the resonance 
is stabilized to be 12.2 eV in our R-matrix scattering calculation. It is unclear why the bump 
in the ICSs of the 5 3 fl„ state is not evident in the previous R-matrix cross sections of Gillan et 



al. ll21ll . Some details of the R-matrix calculations are different in their calculation and ours, e.g., 
they used hybrid orbitals with Slater type functions, whereas we employed SA-CASSCF orbitals 
with Gaussian type functions. These difference may contribute to the difference in magnitude of 
the 2 A g partial cross section. 

In this study, we employed the fixed-nuclei (FN) approximation. As we can see in figure HJ 
equilibrium bond lengths of the excited N 2 states are longer than that of the ground state. Thus, in 
principle, it would be desirable to include the effect of nuclear motion in the R-matrix calculation. 
Use of the FN approximation may be responsible for several discrepancies between our calculation 
and experiments, including bumps in the ICSs of the A 3 E+ and B 3 Tl g states, the position of the peak 
in the ICSs of the C 3 n u state. Although the calculated DCSs agree very well with experimental 
results in general, our DCSs of the A 3 2„, B 3 Tl g , w l A u and C 3 II U states at 12.5 eV are 2-4 times 
larger than experimental results. These deviations in the near- threshold DCSs can also be related to 
the FN approximation. In spite of these discrepancies, good agreements are observed between our 
calculation and experiments in most ICS and DCS cases as we can see in the figures. Agreements 
with the recent experimental results of Khakoo et al.|2|] and Johnson et al. |3|] are especially 
impressive. It is possible to include nuclear motion in the R-matrix formalism through vibrational 
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averaging of T-matrix elements or the non-adiabatic R-matrix method, though application of these 
methods will be a difficult task in the presence of many target electronic states. In the future, we 
plan to perform the R-matrix calculation with these methods including nuclear motion. 
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IV. SUMMARY 



We have investigated electron impact excitations of N 2 molecule using the fixed-bond R-matrix 
method which includes 13 target electronic states, A 3 E+, B 3 Tl g , W 3 A U , 5 ,3 E„, a' 1 ^, a l Tl g , 

w 1 A u , C 3 II M , E 3 I,g, a" 1 ^, c l Il u and c' l l^. These target states are described by CI wave functions 
in the valence CAS space, using SA-CASSCF orbitals. Gaussian type orbitals were used in this 
work, in contrast to the STOs in the previous R-matrix works. We have obtained integral cross 
sections as well as differential cross sections of excitations to the A 3 E+, B 3 H g , W 3 A U , 5' 3 £~, a' x Yr u , 
a l H g , w 1 A u and C 3 n„ states, which have been studied a lot experimentally but not enough theo- 
retically before. In general, good agreements are observed both in the integrated and differential 
cross sections, which is encouraging for further theoretical and experimental studies in this field. 
However, some discrepancies are seen in the integrated cross sections of the A 3 £+ and C 3 n„ states, 
especially around a peak structure. Also, our DCSs do not agree well with the experimental re- 
sults at low impact energy of 12.5 eV, compared to the higher energies of 15 and 17.5 eV. These 
discrepancies may be related to the fixed-nuclei approximation or insufficiency of higher excited 
target states in the R-matrix model. 
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FIG. 1: Potential energy curves of the N2 electronic states. The equilibrium distance of the X 1 ?,* state, R = 
2.068 ao is used in our R-matrix calculations. 
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FIG. 2: Integral excitation cross sections of the A 3 X+ (panel a), S 3 n g (panel b), W 3 A U (panel c) and B' 3 I,~ 
(panel d) states. Our results are shown in thick full lines. For comparison, we include the previous R-matrix 
results of Gillan et al. 12 IN . Schwinger multichannel results of da Costa and Lima 1350, the experimental 
cross sections of Cartwright et al. 11111 . Campbell et al. U3D and Johnson etal.^. 
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FIG. 3: Integral excitation cross sections of the a' x YT u (panel a), a x W s (panel b), w l A u (panel c) and C 3 n„ 
(panel d) states. Our results are shown in thick full lines. In addition to the experimental ICSs in Fig|2j we 
include the results of Ajello and Shemansky J5], Zetner and Trajmar lioj . Mason and Newell |6], Poparic 
et al. [70, Zubek (8D and Zubek and King 
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FIG. 4: Differential cross sections for electron impact excitation from the N2 X 1 !,^ state to the A 3 X+ state. 
Panel (a): electron impact energy of 12.5 eV, (b):15 eV and (c): 17.5 eV. Full line represents our result. For 



comparison, we include the previous theoretical cross sections of Gillan et al.[21], experimental results of 



Khakoo et al. Ol , Brunger and Teubner [ 12], Cartwright et al. Ill ill . Zetner and Trajmar [10] and LeClair 
and Trajmar 13711 . 
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FIG. 5: Differential cross sections for electron impact excitation from the N2 X state to the B 3 U g state. 
Panel (a): electron impact energy of 12.5 eV, (b):15 eV and (c): 17.5 eV. The results of Schwinger multi- 



channel calculation by da Costa and Lima IB5I1 are also shown in the panels. Other details are the same as 
in Figg] 
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FIG. 6: Differential cross sections for electron impact excitation from the N2 X 1 !,* state to the W 3 A U state. 
Panel (a): electron impact energy of 12.5 eV, (b):15 eV and (c): 17.5 eV. Other details are the same as in 
FigfH 
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FIG. 7: Differential cross sections for electron impact excitation from the N2 X state to the B' E M state. 
Panel (a): electron impact energy of 12.5 eV, (b):15 eV and (c): 17.5 eV. Other details are the same as in 
Figffl 
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FIG. 8: Differential cross sections for electron impact excitation from the N2 X 1 !^ state to the a' E~ state. 
Panel (a): electron impact energy of 12.5 eV, (b):15 eV and (c): 17.5 eV. Other details are the same as in 
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FIG. 9: Differential cross sections for electron impact excitation from the N2 X IJ state to the a l H g state. 
Panel (a): electron impact energy of 12.5 eV, (b):15 eV and (c): 17.5 eV. Note that the DCSs are shown in 
logarithmic scale. Other details are the same as in FigQJ 
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FIG. 10: Differential cross sections for electron impact excitation from the N2 X state to the w A u state. 
Panel (a): electron impact energy of 12.5 eV, (b):15 eV and (c): 17.5 eV. Other details are the same as in 
FigfJ 
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FIG. 11: Differential cross sections for electron impact excitation from the N2 X state to the C 3 IT„ state. 
Panel (a): electron impact energy of 12.5 eV, (b):15 eV and (c): 17.5 eV. The experimental DCSs of Zubek 
and King [9] are added. Other details are the same as in Fig]4] 
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TABLE I: Division of the orbital set in each symmetry. 



Symmetry 




B 2u 


B3u 


Big 


B\ u 


B 3g 


Big 


A u 


Valence 


l-3a g 




lb 3u 




l-36i„ 




\b 2g 




Extra virtual 


Aa g 


2b 2 u 




\b lg 




2b 3g 


2b 2g 




Continuum 


5-39a g 


3-35b 2u 


3-35b 3u 


2-\lbu 


5-37b lu 


3-m 3g 


3-m 2g 


2-\la u 
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TABLE II: Comparison of the vertical excitation energies. The present results are shown with the previous 
works of Gillan et al. 12111 . multi reference coupled-cluster (MRCC) results of Ben-Shlomo and Kaldor 13311 
as well as experimental values quoted in Ben-Shlomo and Kaldor 13311 . The unit of energy is eV. 



State 


This work 


Previous R-matrix 


MRCC 


Experimental values 




0.00 


0.00 


0.00 


0.00 




7.89 


7.63 


7.56 


7.75 




8.54 


8.54 


8.05 


8.04 


W 3 A U 


9.38 


9.11 


8.93 


8.88 


a l Ii g 


9.85 


9.89 


9.27 


9.31 




10.06 


9.83 


9.86 


9.67 




10.69 


10.41 


10.09 


9.92 


w l A u 


11.01 


10.74 


10.54 


10.27 


c 3 u u 


11.64 




11.19 


11.19 
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